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We show that it is possible to extract the weak mixing angle a via a measurement of the rate 
for B^ -^ 7r^e"'"e~. The sensitivity to cos a results from the interference between the long and 
short distance contributions. The short distance contribution can be computed, using heavy quark 
symmetry, in terms of semi-leptonic form factors. More importantly, we show that, using Ward 
identities and a short distance operator product expansion, the long distance contribution can be 
calculated without recourse to light cone wave functions when the invariant mass of the lepton pair, 
q^ , is much larger than Aqqq . We find that for g^ > 2 GeV'^ the branching fraction is approximately 
1 X 10^*|Vtd/0.008p. The shape of the differential rate is very sensitive to the value of cosa at small 
values of q^ with dV/dq^ varying up to 50% in the interval — 1 < cosa < 1 at g^ = 2 GeV^. The 
size of the variation depends upon the ratio Kis/Vm. 



Great effort is presently being expended in attempting 
to understand the origin of CP violation. It is hoped that 
in the next generation of experiments we will be able to 
determine all the parameters in the CKM matrix. How- 
ever, the extraction of these parameters is, in general, 
hindered by our inability to perform first principle cal- 
culations of rates due to the non-perturbative nature of 
the long distance QCD effects. A particularly nettlesome 
extraction is that of the angle a in the unitarity trian- 
gle. The standard proposal for the extraction of a from 
B -^ TTTT is hindered by so called penguin pollution, which 
can only be overcome through a cumbersome SU{3) anal- 
ysis. Here we propose to extract this angle via a mea- 
surement of the rate for the rare decayn B — > 7re^e~. It 
is usually assumed that the rate for this process is domi- 
nated by the short distance transition b — > de'^e~ except 
when the invariant e^e^ mass, q^ — {p^- +Pe+)^i is 
of order of charmonium resonances where long distance 
contributions are important (here and below by short dis- 
tance transition we mean contributions to the amplitude 
that are effectively local at distances larger than the elec- 
troweak scale 1/Mw)- However, there is a long distance 
contribution, which arises through weak annihilation di- 
agrams, like the one in Fig. y, which can contribute sig- 
nificantly. The short distance amplitude is proportional 
to Vf^Vfl whereas the long distance annihilation graph is 
proportional to V^j^V*ij. Thus, the interference of these 
contributions leads to a rate which is sensitive to the 
value for cos a where 
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Even a crude measurement of cos a would be of value 
since it would remove a two-fold ambiguity in extractions 
of a from sin 2a. Naively, one would think that any hope 
of extracting a in this way is doomed by the fact that 
long distance contributions are notoriously intractable. 
However, in this paper we show that this weak annihila- 
tion can be calculated in an expansion in l/m(, and as 
when the invariant mass of the electron pair q^ is larger 
than Aqqq. Moreover, as will be seen below, the rate is 
independent of the valence, as well as higher twist, wave- 
functions of the pion, thereby reducing the uncertainty 
in the calculation. 

Let us naively estimate the relative importance of the 
short distance amplitude and of the long distance weak 
annihilation amplitude. The former must involve an elec- 
troweak loop, so it carries a factor of l/167r^, while the 
latter should be suppressed by wave functions at the 
origin, fsf-K/Mg, where fs and Mb are the B- meson 
decay constant and mass, respectively, and /^r is the 
TT-meson decay constant]] In addition, the CKM fac- 
tors are different. So the ratio of the long distance to 
short distance amplitudes is expected to be of the order 
of \Vub/Vtd\l6TT^fBf-K/Ml - 0.07. Here we have used 
/b = 170 MeV, frr = 130 MeV, and \Vut/Vtd\ = 0.6. 
Thus, weak annihilation can easily give a correction of 
10% to the rate, possibly larger. A more quantitative 
calculation is well motivated. 

Let us now undertake a systematic calculation of the 
weak annihilation amplitude in B^ -^ 7r*e+e^, which 
will be the primary focus of this paper. The relevant 
effective Hamiltonian for the long distance, weak annihi- 
lation contribution to the rate for B^ —>■ 7r^e+e^ is 



* It is also possible to use the mode B — > pe"* 
be discussed in a separate publication |1|. 



' ^At small q^ the short distance amplitude is further sup- 

pressed by q'^/Mg. 




FIG. 1. Weak annihilation diagram underlying the decays 
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There are three other diagrams 



with the photon emitted from any of the three light quarks. 
Photon emission from the VF-boson is suppressed by Gf- 
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KbKTd {c{pi/Mw)0 + c'{pi/Mw)0') , (2) 



= uYP-b d-f^P-U 



O' = u-f^P-T^b d-f^P^T^'u, 
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P± = (1 ±75)/2 and T" are the generators of color gauge 
symmetry. The dependence on the renormalization point 
/i of the short distance coefficients c and c' cancels the 
/i-dependence of operators, so matrix elements of the ef- 
fective Hamiltonian are /^-independent. At next to lead- 
ing log order, using Aq^j-, = 225 MeV, the coefficients at 
^ = 5 GeV are @ c = 1.02 and c' = -0.34. 
The B~ — > 7r~e~'"e~, decay rate is given by 
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where (.^ — u{pe-)j'^v{pf,+ ) is the leptons' electro- 
magnetic current, q = pB + Pe- and t = (pu + Pe+)'^ = 
{pb —Pe-Y- A sum over final state lepton helicities is 
implicit. The non-local contribution to the hadronic cur- 
rent h is 
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The two body kinematics with an energetic massless fi- 
nal state hadron is known to factorize [pj in the sense 
that soft gluon exchange between initial and final state 
hadrons is suppressed by l/mh or as{mii). This factor- 
ization results from the fact that the energetic outgo- 
ing light quarks form a small color singlet object which 



does not couple to leading order in the ratio k/Eg ( "color 
transparency" H|), where k is the soft gluon momentum 
and Eq is the energy of the outgoing quarks which scales 
with TTifc provided q^ is not close to g^ax = {Mb — JtItt)^- 
Furthermore, by the same reasoning, the color octet op- 
erator does not contribute to leading order in k/Eq once 
the final state is projected onto the color singlet channel. 
Notice that this is the simplest exa-mple of Bjorkens' color 
transparency argument, since there is only one hadron in 
the final state. There are no assumptions needed regard- 
ing the behavior of wave functions. 

Returning now to our result we find that factorization 
leads to 



+ {7^\jxiOm{0\TijU^)j\0))\B)] (7) 

where k = 2V2G_f Vub^ud^, j\ = dj\P-u and J^ = 
u'y^P-b. 

The first line in Eq. (0) can be computed using an 
isospin Ward identity. The B-momentum acts as a 
derivative on the T-ordered product which then gives the 
matrix element of a commutator. It makes a contribution 
to h^' oi-enf^fePB- 

It is remarkable that the result is not sensitive to the 
internal structure of the pion. The first term in Eq. (M) 
may be written using a lightcone OPE, as a weighted 
integral over the pion valence wavefunction plus higher 
twist corrections M . If one were to perform this calcula- 
tion, one would find that to all orders in as the Wilson 
coefficients are independent of the momentum fraction 
carried by the light quarks. The generality of our re- 
sult allows one to see immediately that in addition it is 
independent of all higher twist wavefunctions as well. 

The second line in Eq. (0) is also insensitive to the 
lightcone wavefunction of the B meson, in that it can be 
computed using a short distance OPE y,g| and a heavy 
quark expansion (HQET) for the 6-quark, provided g^ ^ 
A2 



To leading order we find 



h^ = -feK/^/BP^. 



(8) 



The coefficient of the singlet operator (Q) is a function 
of renormalization point, /i. However, once we have 
matched to the HQET the /i-dependence of the coeffi- 
cient cancels that of the decay constant of the i?-meson 
in the HQET. This invariant combination is the physical 
decay constant /^, or rather, the leading approximation 
to it in an expansion in l/mi,. 

As we will see the sensitivity to cos a is greatest at 
small q^ where the short distance contribution is sup- 
pressed due the rapidly falling form factors. Thus we 
would like to be able to trust our results to as low a value 
of q- as possible. If k denotes the momentum of the u- 
quark in the _B-meson, the OPE is a double expansion in 



q ■ k/q^ ^ AQcunib/q'^ and P/q"^ -^ Aq^^q/q^. The lead- 
ing correction, of order q ■ k/q^ can be computed. The 
remaining corrections are of order 20% for q^ > 3.5GeV^ 
and can be computed in terms of matrix elements of local 
operators. Including the leading correction we find 



h" 



^^UfBpU^ + ^Amt/q 



(9) 



where A = Mb — rnt, is the meson mass in HQET. A 
recent extraction g gives A « 330 MeV. 

We now combine the long and short distance contri- 
butions to the amplitude for B^ — > 7r^e+e^. At leading 
ordcm the short distance contribution to the amplitude 
is obtained from the effective Hamiltonian 110(1 



Atcff = — ^ Vth Vf 
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At leading-log, with as{Mz) — 0.12 and rrit = 175 GeV, 
C7(mb) = 0.33 and C'io(to6) = 5.3 0. At next to 
leading-log Cginib) = —4.3 ||^. There are additional long 
distance contributions from the operators O and C, and 
from the corresponding operators with the u-quark re- 
placed by a c-quark, in which a photon is emitted from 
a u- or c-quark loop. These contributions can be incor- 
porated into a shift in Cg, 



Cg = Cg -I- (c -I- |c')g(TOc/mb, f) 



9(0, f)] (14) 



where q^ = q^ /m^ and the function g is defined in 
Ref. llfl] . The short distance amplitude is given in terms 
of form factors, f± and h, defined by 

{^{p')\drb\B{p)) = {p+p'rf+ + {p-p'Tf- (15) 

{-K{p')\da^"'h\B(j>)) = 2i/(,(pV -pV) (16) 

Including the long distance contribution to the ampli- 
tude, as calculated above, and neglecting the mass of the 
electron, the rate for B — > Tre+e" is 



■''The penguin operators Oa-Og have small coefficients and 
we neglect them. The gluonic magnetic moment operator, 
Os, does not contribute to -B ^ iill at this order. 



FIG. 2. Differential branching fraction, tb dV/dq^ for 
B -^ 7re"'"e~. The solid (blue), dashed (red) and dotted 
(green) lines correspond to cos a = 0,-1,1, respectively. 
The shaded region shows the uncertainty in our calculation. 
For the short distance contribution the BK form factors of 
Ref. tfll have been used. 
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The contributions from the weak operator d'y^bcq^c will 
be poorly described by the function g when q^ corre- 
sponds to the mass of a charmonium state, so we restrict 
our analysis to (f' < rn^. 

We see that the interference effect is largest at smaller 
values of q'^ where the form factors are suppressed. A 
part of the uncertainty in the extraction will be depend 
on our knowledge of the form factors. The form factors 
/-!_ and h may be extracted from measurements of the 
scmi-leptonic B decays and B -^ K*'y, respectively. Al- 
ternatively, we may use heavy quark symmetry and use 
the relation 



Hq') 



(/+(g2)-/_(g^)) 
2mii 



Oil/mb). 



(18) 



Presently, there exists lattice QCD determinations at 
larger values of q^ , where there is little hadronic recoil. 

Given, that presently we do not know the form factors 
away from largest values of q^, for illustration purposes 
we will use the BK model (of Becirevic and Kaidalov 



jni) which satisfies the unitarity sum rules bounds ^ 
and fits lattice determinations at large q^: 



(BK) 



/ 
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N{l'f3) 



(l-g2)(i_^g~2) 



/« ^ ^^^i^, (19) 



1 -7(72 



= 0.54, 7 



g2 



0.8, iV = 



and the 
0.6 and 



parameters are f3 
Mb' ~ 5.325 GeV. Here we have changed the value 
of iV = 0.8 given in Ref. ||l|], which is appropriate for 
TOf, — 2.6 GeV, for the more realistic value of iV = 0.6. 
The value of 7 is fixed by the Callan-Treiman relation 
/o(™b) — fs/f-K- The BK model does not give the form 
factor h. We calculate h using the heavy quark spin sym- 
metry relation (|l8|). 

In Fig. |2| we plot the rate of Eq. ( p7| ) as a fraction 
of the total width T = Tg^ as a function of q^. We 
have used the coefficients c, c' and Cg at next to lead- 
ing order, and the rest at leading log order. We have 
restricted the plot to q^ > 1.5 GeV^ for our approxima- 
tions to remain valid, and to q^ < 8.5 GeV^ to avoid 
contributions from charmonium resonances. For illustra- 
tion we have used \V^i,V*^\ = 0.004, \Vt^V;^\ = 0.008, and 
fg = 0.17 GeV. The solid (blue), dashed (red) and dot- 
ted (green) lines correspond to cos a = 0,-1,1, respec- 
tively and the shaded region represents an uncertainty 



rate and spectrum would constrain the unitarity triangle, 
e.g., a region of the {p, rf) plane Q . The analysis requires 
knowledge of decay form factors. Semi-leptonic decay 
spectra will determine the combination |l^„b/+|. Such a 
measurement can be incorporated in our analysis 0] and 
constrain the unitarity triangle in a meaningful way even 
if separate knowledge of | Vub \ and /+ is lacking. 
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